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SUMMARY 

Studies were made on ae thrust bearing with two 
peer ent Slider boundary conditions, namely, constant 
temperature and adiabatic surface. This study was based on 
the transient momentum and energy equations. A number of 
numerical simulations were performed using different bearing 
ratios and Prandtl and Eckert numbers. 

It was found that the location of maximum temperature 
depended on the bearing ratio. For a bearing ratio of 2.2, 
maximum temperature was near the slider. For a bearing 
ratio of 1.5, maximum temperature was near the stationary 
surface. Temperature computations were carried out until 
transient values were acceptably close to the steady-state 


values. 
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NOMENCLATURE 


a = acceleration of slider 

a = eho, /(Br*)} Sidamenstonlesseacceleration 

Cc aiiewarhgs | = coefficients used in numerical method 

B = length of slider in direction of motion 

Cp = specific heat of lubricant at constant pressure 
E = (Bp /hs) ?7(CpitalescaxtyperdainEckert number 
h = film thickness 

h, = inlet film thickness 

ho = outlet film thickness 

hi = hyhog= dimensionless film thickness 

pga, ok = numerical method equivalents of x, y, and t 


respectively 
K = hy /ho. = bearing ratio 


M,N = number of computational points in the x and y 
directions respectively 


Dp = gauge pressure 

‘ous = pohe/(Bu)* = dimensionless gauge pressure 

P = Cpu/k = Prandtl number 

Q = flow rate per unit width 

en = Qh,/(Bv) = dimensionless flow rate 

ie = time 

oe = pt/h2Z = dimensionless time 

Tix, oy) -Ttemperature of Lubricant eat position Cx y) 
a = constant surface temperature 

TT: = T/T, = dimensionless temperature of lubricant 
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uhé/(By) = dimensionless velocity in 
M=airection 


velocity of slider 

Uhg/(By) = dimensionless velocity of slider 
VELOC. ey “nwy-airec hon 

vho/v = dimensionless velocity in y-direction 
ecordinatve along, the bearing, 01s x (7B 

x/B = dimensionless x-coordinate 

coordinate across the bearing, 0 < y < h 


y/ho = Gimensionless y-coordinate 


thermal conductivity 
absolute viscosity 
kinematic viscosity 


density Ob -lubuicane 
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CHAPTER I 
INTRODUCTION 


Background to the problem 

Numerous studies have been done on _ steady-state 
momentum and energy equations in connection with thrust 
bearings. Hahn and Kettleborough [1] used these equations 
to solve for pressure and temperature distributions ina 
thrust bearing. The effects of inertia terms on pressure 
and load capacity, andthe effect of convective terms on 
temperatures at the surfaces were investigated by Rodkiewic 
and Anwar [2]. Other studies concentrated on varying temp- 
eratures at the surfaces and different bearing ratios [3]; 
keeping the two surfaces adiabatic [4]; a different numeri- 
cal procedure for calculation of flow in finite width bear- 
Pie eo VE MGe ul 5.1.5 

There have been a number of studies done on the tran- 
Sient momentum equation. One of the earliest works was done 
by Ladanyi [6] who investigated the effects of acceleration 
on performance characteristics in slider and journal bear- 


ings. Lyman and Saibel [7] did an -analytical study of the 


N 


transient pressures in an infinite slider but did not assume 


that acceleration was proportional to the distance across 
the bearing as in the previous paper. A more recent paper 
by Venkateswarlu and Rodkiewicz [8] studied the development 
of transient load capacity and drag for three functional 


variations of velocity. 
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Present study 

This study 1S concerned with transient temperatures 
in a thrust bearing. The solution of the transient momentum 
equation is used to solve for the temperature distribution 
in the bearing. As far as this author knows, this 1s prob- 
ably the first study made on the transient energy equation. 
This research was done for a specific range of parameters. 
Therefore, it 1S recommended that further studies be made 


on the transient energy equation. 
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CHAPTER II 


TRANSIENT MOMENTUM AND ENERGY EQUATIONS 


Governing equations 


The thrust bearing configuration considered for this 
analysyS 1S shown in gage; 1. The ‘governing equations, 
assuming an incompressible fluid of constant properties, 


or the two-dimensional bearing are 


pom toc + vou ie op: tind oe (1) 
ot Ox dy Ox ay? 
du + dv = 0 Gay 
Ox ay 
pGo ome uo T Sv Te) = Bad OT ous 9) (3) 
ot OX Oy ay? dy 


with the following velocity and pressure boundary conditions 
Cigar ty) em Uke yan 
Wee list) as AX) tae oe ery 0) =n 
(4) 
mie te) B= eB at) = 0 
mepxeONte b= al) t= at 


and the following two cases of temperature boundary 


conditions 
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Fig. 1. Coordinate system of thrust bearing 
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The transient momentum equation (1) with the boundary 
conditions (4) have been used in ref.[8,9]. In these papers 
and in the present study, U = at where a is the accelera- 
tion of the’slider.’ This’méans that the slider is station- 
ary when the time is zero. This study will be concentrated 
on the transient energy equation (3). 


Equations (1) and (2) can be combined and together 


with the energy equation (3) can be written in dimensionless 


form as follows 
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The corresponding boundary conditions become 
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quantity h*(x*) describes the shape of the stationary part 
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of the thrust bearing with reference to the moving slider. 


The shape assumed is for a linear bearing, namely 


RicGxesla = Ket: (1 aK) xe (12) 


The solution of equation (7) yields the velocity 
GrSErIbDUCHON Whawcthe bearindgatjtime t”.. Since the fluid 
properties are constant, the solution of equation (7) is 
uncoupled from equation (8) and can be solved independently. 
The VelociEey aistributiony Sak time it’ 91s used im equation 
(8) bo Solve for the temperature distribution in, the oil 
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Numerical method 


The dimensionless differentials in equations (7) and 


Caieean be written in fLinivte-difterence torm as*fclilows 
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respectively, and x* has been replaced by i, y* by j, and t* 
by k.‘  Theforigin is fixed with respect to the stationary 
surface. yw Bmitially, .u°(i),53 Ab) -= 0 for the entire bearing. 
wher k => the bearing acquires a velocity u*(i;3,2) eat" 
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Fig. 2. Computational molecule (x position is fixed) 


k = 3 from equation (13). In the next step, the dimension- 
less flow rate Q, corresponding to 
ht 
og =f uray’ (16) 
0 
1s obtained using Simpson's rule on the values u*(M,1,k), 
oy TOe(MAN [ka neThema~tieis decreased by’ 1, 1 =(1-- 1,°and 
EhemoOOMeUtatTOn OF ti Cimels hw prea. bu GueN ike cand or is 
performed assuming a value for (dp*/dx*),. It is known that 
GO teat teorrect,  \shouldsbeaetunction of /Lime. only and must 
NOt aaveWwieh eChalnigiind, peor sa . Therefore, the computed 
flow rate at 1 and M must be equal. If there 15 a signifi- 
cant difference, a new value of (dp*/dax*); is used. The 
pressure gradient at 1 1S adjusted until a satisfactory 
Q7 is obtained. It is found that the relationship between 
dp*/dx* and Q* is approximately linear. Therefore, a better 
approximation of dp*/ax* can be obtained. The final value 
of pressure gradient at i is found by iteratively correct- 
ing the previous value. 
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0 
is per formed using Simpson's rule over the M - 1 intervals 
of 1/(M- 1) each. It must be equal to zero in order to 
Sattsty the boundary conditions (9). Lf sibs snotnso,, tne 


value of (0p*/dx*)m must be adjusted and the computation 
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of u*, Q* and dp*/dx* is repeated. This process is repeated 
iments the wniLeqral wsteloseltoezero.» It is also found)that 
the relation between (dp*/dx*)m and pf iS approximately 
linear and therefore a better approximation for (dp*/dx*),» 
is possible. The pressure distribution can be determined 


from the following equation 


+ 


Danie ko) an CMe a elriek) Vk (18) 
Ox 
0 

Once the velocity distribution is known at k = 3, the 
temperature distribution can be calculated from equation 
ia veh its coeriicients. in equation (15). 

The temperature distribution is solved iteratively by 
using the Gauss-Seidel method with relaxation. Equation 
(14) together with particular Prandtl and Eckert numbers 
are solved by this method. Iteration 1s stopped when the 
present set of values are in close agreement with the 
previous set. The calculations are based on the assumption 


that one upstream solution is known. Formulation of 


equation (14) is based on a 4-point computational molecule 
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Therefore, the calculations proceed from the entrance 
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Whetve l=) )etoO ithe exit where 1°=M, in -one unit increments. 


TnemeLcotmocre Ol Tol, we Por eye= |, gee, N? iS). known, 
as shown in. equation (10). Note that the computation of 
temperature goes from the entrance (i = 1) to the exit 


(i = M) which is exactly opposite to the way velocity is 
Calculated, After a satisfactory set of temperatures have 
been calculated, k can be incremented and the computation 
can be repeated for the new time interval. 

This procedure has been programmed in double preci- 
Sion on the Amdahl 470/V8 computer. The simulation has been 
performed for a number of bearing ratios and Prandtl and 
Eckert mumbers. A single simulation of the preceding 
problem took approximately 130 seconds of central process- 
iidubimes s)he lubrication sol! wsed has tthe following prop- 


erties (at 40°C): 


Cpe et. COaki/ko K = (O01 TS Wm 2k 

i= 0.000286 Ns/m-* (=) O45 en Oem eS 

op =8O7a ko/me p= 119 A 
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The dimensions of the thrust bearing are as follows: 
JB = 0.0914 m ne o462 x 10 om 


The Eckert number can be calculated from the above data. 
A comparison between the type of Eckert number used here 


and other Eckert numbers used in ref.[1,2,3] would show 
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that the slider velocity, U, does not appear in the former. 
For this study and others [8,9], variables that are used 
in dimensionless constants must not ‘vary ‘with “tame, 
Hence, the slider velocity which varies with time has been 
replaced by Suitable time-invariant variables. The comb- 


imations used are shownrin-Table I. 


Constant slider temperature 
Prandtl Eckert 


1S) 0.08140874 


be, 0.08140874 


Adiabatic slider 


Prandtl Eckert 
0.08140874 


0.08140874 


Tablee!lli Bearingetatsesm BPrandti and Eckert 
numbers used in simulations 
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CHAPTER III 


STEADY-STATE MOMENTUM AND ENERGY EQUATIONS 


Governing equations 


The thrust bearing configuration considered in the 
Steady-state analysis of momentum and energy equations is 
Similar to Fig. 1 except for the absence of the time axis. 
The governing equations for a two-dimensional bearing of 


constant properties are as follows 


pueliow Vous) p= = sdprripcru (49) 
OX Oy Gx ay" 
Gu OV e= 80 C20) 
Ox dy 
pepieuol + wel) Meo Fe enewour) = G20) 
OX dy Oy? y 


with the following velocity and pressure boundary conditions 


chery ee ys) kaa Cie Or ae 
G22) 

ws, Oye 10, pCO) = p(B) = 0 
ana the ~following two cases ‘of temperature boundary 
conditions 

aor ee Doh) Fe Op ae eal 23) 
or 

T(0,y) = T(x,h) = 7,, aT = 0 (24) 

ay ity = eo 


Equations (19) and (20) can be combined and together 
with the energy equation (21) can be written in dimension- 


less form as follows 
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(2%) 
iG eel) = All Oeeoeseo. (hiar= 0 
and 
Oey ety ec (RIS eee ec be) pic (28) 
Ou 
mCGey oT. (xs ih ee ay 0 = 0 (29) 
dy" ly*= 0 


ime “the, Sboundary condita ons’ (27), (28), (29), ‘the 
Quantity h*(x*) describes the shape of the stationary part 
of the thrust bearing with reference to the moving slider. 
The shape assumed is for a linear bearing aS given by 
equation (12). The solution of equation (25) is uncoupled 
from equation (26) because of the assumption that fluid 
properties are constant. Therefore, equation (25) can _ be 
solved independently to yield the velocity distribution in 
the thrust bearing. The results are used in equation (26) 


to solve for the temperature field in the lubricating film. 
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Numerical method 


The differentials in equations (22) and (23) can be 


Wereeenm@in ornite differencesrorm as follows 
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The computational grid is the same as before. For 
the steady-state calculations, the time variable does not 
appear. An initial set of starting values for the velocities 
in the bearing was generated from a program using equation 
(22) without the inertia terms. These values will serve as 
the initial guess to the full equation (22). To start the 


iterative process, an arbitrary value of (dp*/dx*), 1s 
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guessed. The value used must be negative because of the 
boundary condition (22). 

Using the initial pressure gradient at M, and Sstart- 
EnoawtthneMperM rine (M, 7). 9)=) 1, 24.., N the initial velocity 
estimates can be calculated from equation (30) using the 
Gauss-Seidel method with relaxation. The iteration 1s stop- 
ped when the present set of values are in close agreement 
With the previous set of values. The rest of the numerical 
procedure is quite similar to that used for the transient 
equations. Any references to time is ignored. 

The final velocity distributions are saved on compu- 
ter files. These values are then used in the energy equa- 
tion (31) together with the appropriate Prandtl and Eckert 
number (Table I). The iterative technique used 1S again 
the Gauss-Seidel method with relaxation. Iteration is 
terminated when the present set of temperatures are in 
close agreement with the previous set. 

A single simulation of the steady-state equations 
took approximately 20 seconds of central processing time. 
The programs used to solve the equations have been written 


in double precision on the Amdahl 470/V8 computer. 
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CHAPTER IV 
DISCUSSION OF RESULTS 


Transient temperature profiles for the case of constant 
slider temperature . 

Grapns have been drawn for bearing ratios ,K, of 14.5 
Soe. eT aISCUSSION wubl) “apply “co. both oratios. “If 
there are any Significant differences between the two 
ratios, the discussion will be separately applied to! each 
ratio. In the following discussion, the reference to time 
1S Same aS Saying that it 1S dimensionless time, t*. 

The slider is uniformly accelerated from rest when 
the time 1S zero. When a time of 1.0 1s reached, the slider 
assumes a uniform velocity of 1.0. During the acceleration 
of the slider, there is a sharp temperature rise close to 
the moving slider surface. This is clearly seen in the 
Die ee PhOn =) 0.4, “O.e) Were ‘times: st] 0.5 sand "150 
CO See age a |S?) The high temperatures near the 
moving surface are due to viscous effects. At a time of 
t*= 0.5, viscous effects have weakly penetrated across 
(“Ge lfrom “the “slider to ‘the ‘stationary surface } the 
beecing. sHowever, ithe resulrani velocity proriles as seen 
(ish LCSteisee los rand 815. sehow het thel velocity: oF 
fluid closer to the slider is more fully developed than 
fluid “adjacent © to. the stationary surface (‘or pad ). 


Therefore, the highest shear stress occurs at the slider 
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Surface. This is the case as shown in the plots of shear 
Stress across the bearing (Srigs.e |Cisaae 18% 238 — pao et 
Thus, the highest temperatures will be located near the 
foyimg-sesmriacez Atuthe nexechtame tintervadl? (care an fhe 
vescomsmetiects haverystromgly “penetrated tthroughoutemthe 
bearing. This iS easily verified by observing that the 
velocity profiles at this time are not too different from 
then@steady=statelsprofaiese “(eFigs, 148 16. )« The temp- 
erature profiles still exhibit sharp temperature gradients 
evesetiiottheasliderghtThis fPasitdue toeSthe facte that the 
slider is Still being uniformly accelerated. Fluid 
particles in the bearing do not have a chance to catch-up 
to the steady-state conditions. The temperature profiles 
atetherexite{ox*=51.0 Jt iaresquite different fromithe. other 
upstream profiles. This 1S due to the effect of pressure 
on the development of velocities in that region. The 
pressure increases from zero gauge pressure at the entrance 
ameereaches alimaxamumi@in hes bearings*sAftertbthate pointy 
the pressure decreases back to zero gauge pressure at the 
exit. Pressure distributions are shown in Figs. 27 and 28. 
Velociiyiprofidessiocatedsslightiysafter the pointeofry zero 
pressure gradient and earlier are adversely affected by the 
pressures present. This causes the velocity profile to 
assumettavebulcey gnebthes directionpect wilow.f Attay point 
approximately 0.1 units beyond the zero pressure gradient, 
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An analysis of the full momentum equation (1) will give the 
E@aSonmbopestatingtthabethertvelocitytisvbinéar sat iacrpoint 
beyond the zero pressure gradient. If it is assumed that 
Steady-state is almost reached, the first term on the 
left-hand-side is very small compared to the other terms. 
Inertia terms will remain on the left-hand-side of the 
equals sign and pressure gradient and viscous terms on the 
right-hand-side. In order to get a linear velocity profile, 
the pressure gradient must balance the inertia terms. 
Therefore, a Slightly negative pressure gradient 15 
necesSary to eliminate the inertia terms. When the inertia 
terms are neglected, the case of zero pressure gradient 
gives a linear velocity profile. Beyond that location, the 
pressure favourably affects the velocity profiles. They 
acquire a bulge in the direction of flow. Shear stresses 
aeestnitiallyshighest«dt the moving suriace® tUEIgs+. '49n-a27 
23, 24 ). When the velocity is linear, shear stress is the 
Same across the bearing. This can be seen in the velocity 
(AGUOe= He) Randishedrostress MOFigqye25the distributiéonstaat 
ties Uey Beoreba Ibearangeratio oped Ste Rhiterutha tapos: tion} 
the highest shear stress drifts to the stationary pad 
( Figs. eo) 23 )emorThereforeyeheatSgqenerationiwitly bet tthe 
highest at the stationary surface. However, the bulge in 
the exit velocity profiles do not appear until t*= 1.0. 
Whenatheatime rebsds5, 2the Gsheartstrésstatethenexut 
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top half of the bearing. Hence, heat generation is about 
the same across the upper half of the bearing but increases 
in the lower half and reaches a maximum at the slider. 
There is a slight bulge in the temperature profile near the 
Moving agsunfaces «EF irgs¢K6,640 7): The profile gradually 
tapers off to the constant temperature boundary condition 
atiyehe Stationary pad. 6 lAmethen4aneut fatime of «1 20;,a second 
bulge in the temperature profile can be seen near the 
Stationary pad. There 1s a definite shift of maximum shear 
Stress from the slider to the stationary surface ( Fig 22, 
26@)aretThars «shifthis néeflectedWin the temperature plotsscf 
both bearing ratios. The presence of two local temperature 
maxima at this time could be explained in the following 
manner. Initially, shear stresses, from the entrance to 
the point of linear velocity, are the highest on the moving 
Sugrace ROEM GS Ss 1+Ftm-ie2heu2se-425e) sheASdal result, the chigh- 
est temperatures are developed near that surface. When the 
maximum shear stress shifts to the stationary surface, the 
highest heat generation is now located on this surface 
( Figs. 22; 26). Therefore, local temperature peaks start 
to form near the stationary Sn ee Ln sother (words, the 
development of temperature peaks at the stationary surface 
occurred in the latter part of the bearing. When the time 
is 1.5, the two ratios are discussed individually. 

For apbearingoratio.wof 2.12; the ssheamrstresses yhni ithe 


bearing correspond closely to the location of maximum 
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temperatures. The temperature profiles at x*= 0.1, 0.4 and 
0.7 have maximums which correspond well with the location 
SsGeiaxr@mumrmshSar Shressa@y( Figs. S31605,\1954u270) % The 
highest shear stress at these locations are at the moving 
Surface. At the exit, the correspondence does not hold for 
EHES PaLLtOR( “Pigs. 76, *22t).Viinsteadys*thethighest etemperas 
tures are near the slider and maximum shear stress is at 
the stationary pad. The higher shear stresses at the 
pad result in higher heat generation but the overall 
temperature increase 1S not great enough to exceed the 
temperatures at the slider. When a comparison is made 
beeween the@qraphs Sat" xtae ON9rand WOR CRY gswes,ec) yeoit 
can be seen that the temperature rise near the pad is 0.04 
OGri2. 5°64) Sand tat! theermovangegsurfacents 2000164015 .08EN)% 
However, the temperatures near the slider still remain 
higher than the temperatures near the pad. 

Por a) bearingyraticororvi25, the locatron foimemaximim 
Shear stress roughly corresponds to the location of maximum 
temperatures. At x*= 0.1 and 0.4, the temperatures across 
GO%™ of" thePGbearing are about )the same? (GFIqs.°.7)e8r)e but 
the highest shear stresses are at the slider ( Figs. 23, 
24 ), The shear stress at x*= 0.7 is approximately the 
Samenracrossiuethewebearangiy (Figs 26t). “This. astexpectied 
becauser theylocatione ot “Pinéarl ’velocity "(Mrigqavu7e)ciis 
abourro?1tunitsY from Sthe “point%ofe zero’ pressure®* gradient 


( at 0.6 ) as seen in the pressure distribution graph 
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( Fig. 28 ). The temperature at this point is a maximum 
near the pad ( Fig. 9 ). At the exit, heat generation is 
hmogherhar Bheripadiathanayatesthenostidere( Baigen26¢)valThe 
temperatunce ancreasetfromexist0oh# From 120 asrQ@.04 ¢42.57?6) 
Nhearytherpadhwandns0.008 FleaSa6SC )}anear thefisSledersulat as 
interesting to note that the increases in this region are 
about the same for both ratios. However, the temperature 
increase near the slider is slightly higher for K = 1.5. 
When a comparison iS made between corresponding shear 
Stress profiles for the two bearing ratios ( compare 
Brash 19 so22 with Bagse 23aci26 );FiteissapparentSthas the 
overall shear stress in the bearing with the smaller ratio 
is higher than the larger ratio. This means that the heat 
generated in the former bearing is higher and will result 
in higher temperatures. This 1s the case when temperature 
profiles at the same location are compared ( compare 
FagS.a3e7n6liwithrFigs.*7as,101)2 etThe. highestutemperatunes 
fesmbothtbearing rabhoseeceug at the rexitel Figs. GewhOe). 
Een tampratioreofie222, i jtheehigqhestrtempertattreteished.0865 
(65.5°Gs)mand,efer asratio ofpi.5,ithe highéstetemperature 
is 1.12 ( 77.6°C ). The surface temperatures were assumed 
to be at 40°C. Therefore, the maximum temperature rise for 
theyhigher and lower fbearninghrathostarepe25e5°Ceand 37.6°C, 
respectively. According to ref.[11,12], a temperature rise 
of approximately 30°C should not be exceeded. However, any 


substantial deviation from this value will affect bearing 
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performance. Therefore, the higher ratio bearing has an 
acceptable temperature rise. The temperature rise in the 
the other bearing is about 25% above the accepted value. 
This is not too great a deviation from the norm. For the 
lower ratio bearing, the flatness in the profiles could be 
due to a combination of high heat generation and the result- 
ing increase in heat conduction across the bearing. The 
earlier shift in maximum temperature, from the slider to 
the pad, for the bearing ratio of 1.5 1s due to the change 
in favourable pressure gradient at x*= 0.6 instead of 
iano See Lorerthee othercitrakitca shebroserSpeo el aerSiecesthe 
temperatures in the bearing were already quite evenly 
distributed, a change in maximum shear stress is quickly 
reflected in the temperature profiles. 

With reference to the velocity and shear stress, it is 
seen that the profiles from a time of 1.5 to steady-state 
foumnvanasingte Curvyer(thigqsvPiis;se@260). cht! would seemethat 
steady-state has been reached at a time of 1.5. However, 
the temperature fields indicate that steady-state has not 
beecna@teecachedutevenpuattianstime sofon3e0 (i hbgqser3aq 10u), 
Rewoanaon there are differences in ae velocity curves for 
times greater than 1.5 but the differences between the 
values are not significant enough to be evident in _ the 
Gbapist, Thrse¢ustalseithe ecaselforntheushean strnessacurvess 
But these ae differences are greatly amplified in the 


temperature graphs. The largest differences occur at the 
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exit because smaller differences upstream tend to have an 
additive effect on subsequent downstream temperatures. 
Therefore, it can be surmised that transient velocities and 
Shear stresses reach the steady-state condition faster than 
transient temperatures. There 1S a greater temperature 
difference in the thrust bearing with ratio of 1.5 than the 
one with a 2.2 ratio. The differences between transient 
( t*= 350. )sand steady-state temperatures are less than 5%. 
This 1S generally the case for most of the profiles except 
for the profiles at the exit where the difference is less 
than 10%. For the velocity and shear stress distributions, 
the differences between the transient and steady-state 
values#eiseounsignificantuatter .t Hove Se 

Therefore, it can be concluded that steady-state cond- 
Hoaons@exaisturfor bthelivellocitye andi sheatastresstdistri- 
butions when atime of 1.5 is reached. This assumption is 
not valid for temperature profiles at that time. However, 
the profiles at atime of 3.0 adequately represent steady- 
State conditions. The differences could be further reduced 
by carrying out computations for a longer time period but, 
the changes between succeeding temperature values decreases 
ase time ancreasess: -Thereforepe there owill?.not be much 
improvement for computations beyond a time of Soe 
However, the costs. will still continue to increase at a 
constant rate. The law of diminishing returns made it 


necessary to terminate computations at atime of 3.0. 
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Transient temperature profiles for the case of an adiabatic 
slider 
A comparison between the temperature distributions 

for an adiabatic slider and a constant temperature slider 
Shows that the only significant differences occur close to 
thensimderhe Faeosehs < WOgande29--936)) .entheicunsves#fotiow 
very Similar paths at distances greater than 0.3 units away 
from the slider. The high temperatures at the slider 
Susface are the resubt Ofy noa heat inconduectionsethrotughntie 
Slider. Therefore, any heat generated at that surface can 
only be conducted into the fluid. When a surface is adia- 
batic, it 1S usually the case that the highest temperature 
will be found on that surface. However, an examination of 
the graphsewilicrevealsthatythe maximum temperatures in the 
bearing occur in the fluid. The shift in maximum tempera- 
ture from the adiabatic surface could be due to a combina- 
SHON OGbmftactorSmenFacthorsfsuch gasonconvectronmreonduct non, 
ancuvuscousheffectsoncoutdereontni butecatohesuche ladeshr it: 
Perhaps, an analysis of individual energy terms may throw 
some light on this phenomenum. It may also be possible to 
uncover the reason why there are two temperature maxima at 
thetexprrofvithegbearang wrthathe rat resof 4202. 

The first energy terms to be analysed will be the two 
convective terms. The first term 1s concerned with convec- 
tion along the bearing and the other is concerned with 
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will be analyzed ( Fig. 37 ). It shows that there is very 
high convection close to the slider surface compared to 
convection near the pad. This high rate persists over most 
Sor@thenbearingiiength?.(,xt=10.7e)m S8his nwouldycmeantthat 
higher temperatures would be near the moving surface. As 
seen in the graphs ( Figs. 29 - 31 ), the higher tempera- 
tures are located in the lower half of the bearing. Around 
xO= OS OF convectionne irs Saboutithetsame’ -acrosstdeuncinthe 
bearing and is lower for the remaining 30% near the pad. 
Therefore, the temperatures slightly above the middle of 
eheebearingwiivige= 095 Penshouldigbest stark ingnvtoti catch-up 
with temperatures near the slider. This increase is part- 
lally reflected in the temperature curves at x*= 0.7. At 
the exit, convection is very high between 0.5 and 0.7. The 
Eerie eacheson avavaluevotieaboutur0sd tunatsavel fethestemp- 
Fature. proraiesdatt thetexitraceexaminedatiFag. S82 Xhe aitius 
seen that the location of a second temperature peak also 
CcGUrSein whatuiregqion.epfheretore, Gitecanebe concluded that 


the first term does help explain the overall shapes of the 


temperature profiles. The overall trend in the convection 
term has been observed in ref.[2]. In that paper, a plot 
of temperature gradient versus x* was done. The general 


Shape of the curves at the pad and slider were about the 
same. it) is0 stall mot @owitteedcleare whyghthe ‘haghest 


temperature is not located on the adiabatic surface. 
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Now, tethe. other ticonvectionotermi(6Fig. 388) wil be 


discussed. For most of the length of the bearing, this term 


does not contribute much to the overall temperature picture. 


From the entrance to 0.7, this term does not exceed 0.025 
units. inltwisebeyond the pointaof linear lvelocitye(s7=*0rs) 
that this term begins to contribute significantly to the 
overall temperature profile. Convection is higher at the 
top 30% of the bearing. Such high convection rates help to 
rapidly increase temperatures near the pad. This observa- 
tion was made earlier for the case of constant slider temp- 
erature. However, the contribution of both convection terms 
do not make the temperatures near the pad higher than 
temperatures near the slider. 

Therloraph of thesconduction term ( Fig. 39 ) shows 
a somewhat confusing jumble of curves that have been plott- 
ed for various distances across the bearing. At the exit, 
it would seem that the highest conduction occurs around 0.4 
anasto. CONnauetton tapersmomige to, either’ side am) this 
region. In that region, there is the presence of a saddle. 
To either side of the saddle, the temperatures are higher. 
This is.seen in the temperature plot ( Fig. 32 ). 

An analysis of the viscous term ( Fig. 40 ) shows 
that viscous effects across the bearing stay quite constant 
Brom x) =) OletoO, O07) AS »-expectea,) ithe. highest viscous 
effects are close to the moving surface. Beyond the point 


of linear velocity ( about 0.8 ), there is a reversal of 
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viscous effects. The highest viscous effects are now at 


the pad. At the exit, the@iverssiviogh viscous effects (0.22) 


at the pad contributed a great deal 
increase in the region. 

In summary, it can be said that 
contributed in some way to the overall 
The lower temperature on the adiabatic 
to numerical inaccuracies. Numerical 
ture at and near the moving surface 


acclirate fanite difference formula. 


to the temperature 


all the energy terms 


temperature profiles. 


Surface could be due 
Solution of tempera- 


is based on a less 
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CHAPTER V 
CONCLUSIONS 


Studies were made on a thrust bearing with two 
different slider boundary conditions, namely, constant 
temperature and adiabatic surface. A number of numerical 
Simulations were performed using different bearing ratios 
and Prandtl and Eckert numbers. 

For the case of a constant slider temperature, it was 
found that the location of highest temperatures depended on 
the bearing ratio. The highest temperature was located 
near the slider for a bearing ratio of 2.2 and near the pad 
for Wa?e bearing pwatro of 1.5. Steadysstete’ wast) reached 
(withingl0%)yat a time oh G0. 

For the case of an adiabatic slider, the temperatures 


on the slider were much higher than in the previously con- 


Sidered case. However, the overall temperature profiles, 
except close to the adiabatic surface, are not too differ- 
ent from the other case. If the-.temperature at the surface 


is too high, cooling as required. Hence, this case would 


indicate whether a constant temperature slider is needed. 
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Fig. 3. Transient temperature distribution at x*= 0.1 for 
K = 2.2 (constant slider temperature) 
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Fig. 4. Transient temperature distribution at x*= 0.4 for 
K = 2.2 (constant slider temperature) 
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Fig. 32. Transient temperature distribution at x*= 1.0 for 
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